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ABSTRACT 

CORRELATION  MEASUREMENTS  IK  THE 

awHUAM   jffifim.   sS.  „  ¥  MB  f  ilfPiir * 

by 
WILLIAM  L.  MARTIN,  III,  LIEUTENANT,  U.  S.  NAVY 
ARTHUR  T.  WHITE,  LIEUTENANT  COMMANDS* ,  U.  S.  MAW 


Submitted  to  the  Department  of  Nuclear  Engineering 
on  July  21,  1958,  in  partial  fulfillment  of  the  require- 
its  for  the  degree  of  Maater  of  Science, 


The  objective  of  this  thesis  ens  to  assemble  an  instrumentation 
syetest  and  make  United  measurements  of  correlation  quantities  in  the 
boundary  layer  of  a  flat  plate  with  an  unheated  starting  length. 

A  system  consisting  of  suitable  probes,  amplifier,  calibration  and 
measuring  circuits,  and  voltmeter  una  assembled  and  tested,   A  series 
of  measurements  of  the  characteristic  quantities  of  turbulent  fluid 
flow  were  made  by  a  traverse  through  the  boundary  layer  at  one  location 
on  the  plate. 

Although  there  is  considerable  scatter  in  the  results,  the  system 
can  be  expected  to  give  satisfactory  results  after  the  recommended 
equipment  modifications  and  test  procedures  are   effected. 


Thesis  Supervisor:  R.  J.  Nicker son 

Title:  Assistant  Professor  of  Mechanical  Engineering 
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I.  IETKOPUCTIOW 

In  the  theory  of  turbulent  fluid  flow  end  heat  transfer,  the  correlated 
quantities  uv  and  vt  have  particular  significance.  The  simplified  equations 
that  describe  the  boundary  layer  over  a  heated  flat  plate  aay  be  written  est 
(See  Appendix  B  for  derivation) 
Momentum  transfer: 

JL<P+pU       +p»2)  *    |y  </||  "puf    -puv*)j  (1) 

Energy  transfer! 

~<ijT  +  ut>»    §^u||    -W    -vt).  (2) 

Analytical  solution  of  these  equations  is  possible  only  under  very 
restrictive  conditions.  For  such  solutions  the  following  definitions  of 
the  eddy  shear  stress  coefficient,  £a  and  the  eddy  energy  transfer  coefficient, 
fih  sre  Introduced: 

•»  mm 

uv  ■  -€m  ~y  (3) 

and  s 

Assumed  analogies  between  these  coefficients  together  with  suitable  boundary 
conditions  has  been  an  accepted  method  for  solving  problems  for  a  great  many 
years.  Very  little  work  has  been  done,  however,  to  actually  measure  uv  and 
vt  and  experimentally  substantiate  the  theory. 

the  literature  contains  little  data  ether  than  that  of  Johnson  in 
ref.  1.  Be  obtained  measurements  of  these  and  other  quantities  in  the 
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boundary  layer  over  a  flat  plata  with  an  unhaatad  a Carting  length.  Bis 
results  ara  given  for  only  one  vertical  exploration  of  the  boundary  layer, 
mmd   they  agree  with  the  expected  behavior  of  uv  and  vt  as  functions  of 
distance  from  the  heated  plate. 

It  has  been  proposed  that  a  thorough  investigation  of  the  thermal 
boundary  layer  properties  be  undertaken  at  H.I.T.   Therefore,  it  is  of 
fundamental  importance  to  design  an  instrumentation  system  of  sufficient 
sensitivity  and  accuracy  to  evaluate  such  properties. 

It  mas  the  purpose  of  this  thesis  to  assemble  such  an  instrumentation 
system  and  make  preliminary  measurements  over  a  heated  flat  plate.  Theoreti- 
cally a  hot-wire  anemoneter  system  will  provide  sufficient  data  for  calculating 
uv  and  vt.   A  basic  system  was  assembled  and  tested.  Certain  components  were 
available  in  the  laboratory,  and  the  remaining  items  were  selected  from 
commercially  available  instruments.  Appendix  f>  contains  a  description  of 
the  equipment. 
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II.  PROCEDURES 

The  experiment  was  conducted  In  the  Heat  Measurement  Laboratory  using 
one  of  the  small   scale  wind  tunnels.  The  boundary  layer  was  developed 
over  a  40"  aluminum  plate  resting  on  the  bottom  of  the  tunnel.  The  last 
twelve-Inch  section  of  the  plate  was  heated  by  a  coil  of  resistance 
heating  wire  on  the  underside.  Measurements  were  wade  along  a  vertical 
traverse  through  the  boundary  layer  ten  inches  from  the  start  of  the 
heated  section. 

The  hot-wire  anemometer  equipment  consisted  of  hot-wire  probes,  a 
Flow  Corporation  HWt  amplifier  and  bridge  unit  and  a  Ballantlne  320  true 
root -mean- square  volt-meter. 

Two  different  hot-wire  anemometer  probes  were  used.  One  was  a  single* 
wire;  the  other  a  V-wire  probe.  Appendix  D  describes  the  equipment  in 
greater  detail. 

Because  the  V-probe  soon  proved  unacceptable,  the  single-wire  probe 

was  used  thereafter.  Orienting  the  wire  first  at  plus  45  ,  then  at  minus 

o 
45  to  the  free  stream  velocity  gave  the  desired  data.  Since  only  the 

time-averaged  quantities  enter  the  calculations,  data  for  each  orientation 

did  not  need  to  be  taken  at  the  same  time. 

After  thermal  aging  for  a  half-hour,  the  wire  was  carefully  calibrated 

by  recording  the  wire  current  at  each  resistance  ratio  for  various  wind 

velocities.  A  mlcromanometer  was  used  to  obtain  accurate  velocities. 

2    / — 
A  plot  of  I  ws  V  0  appears  in  Fig.  8  for  Jf  ■  1.2  and  1.4.  These 
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curves  agree  with  the  characteristic  straight  lines  predicted  by  the  King 
equation. 

From  these  curves  it  is  possible  to  calculate  1*  from  equation  (36), 
where 

c2 

*■  s  Try* —    •  (5) 

As  equation  (72)  predicts t  the  quantities  uv  and  vt  are  related  by 
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where  the  subscripts  a  and  b  denote  the  quantities  for  the  two  wire  orienta- 
tions.  It  i3  important  that  the  wires  of  whatever  probe  is  used  be  correctly 
angled  to  the  flow,  i.e.*"  45  for  this  case. 

The  procedure  for  measuring  and  calculating  these  quantities  follows: 

Kef.  12,  the  HUB  Instruction  Book,  describes  how  to  compute  the  square- 
wave  current  and  how  to  set  the  desired  resistance  ratios.  Briefly,  the 
"Bridge  Null"  dial  of  the  HWB  is  adjusted  to  balance  the  bridge  for  a 
"cold"  current  in  the  wire  of  about  one  and  a  half  milliamps. 

At  one  of  the  available  resistance  ratios,  the  current  is  increased  in 
the  wire  until  the  bridge  is  again  balanced.  The  values  of  resistance  ratio 
and  "bridge  hull"  dial  setting  then  enter  the  calculation  of  square-wave 
current  as  shown  in  Appendix  E. 

The  following  sequence  of  readings  was  obtained  at  the  maximum  tunnel 
velocity  of  26  ft/sec  with  the  plate  first  unheated,  then  heated: 
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With  the  wire  at  plus  45°,  values  of  M  ,  Mtcn  and  la  were  recorded 
at  resistance  ratios  of  1.2  and  1.4.  This  was  repeated  at  predetermined 
distances,  y,  above  ti»«  plate  as  Measured  by  the  micro-traverse  gear, 
starting  where  the  probe  shank  touched  the  plate. 

Values  of  %n,  !y-acn  an^  lb  we   subsequently  recorded  with  the  wire 
at  minus  45  for  the  same  values  of  y. 

The  purpose  of  the  hot  and  cold  plate  readings  becomes  clear  when  one 
assumes  t  =  0  for  the  unheated  plate.  Then  equation  (72)  is 


uv 


=  -if!*) 


(7) 


'sea  -  ''sn  b      ^cn  "  "sn  *  g 
However,  with  the  heated  plate  the  full  equation  applies. 

The  remaining  part  of  the  coefficient  of  vt.   ■?■■  r jy  "—  ■  '•■ 

may  be  calculated  at;  shown  In  Appendix  . 

leaking  ocusuremats  at  tvro  resistance  ratios  gives  t**o  simultaneous 
equations  at  each  y.  Solution  of  these  equations  gives  uv  and  vt  for  the 
hot  plate,  and  uv  for  the  cole  plate.  With  the  plate  cold  uv  is  over- 
detenainec,  so  the  average  of  the  two  equations  was  t^ken. 

A  traverse  of  the  boundary  layer  with  the  wire  perpen  dlculat  to  the 
swam  flow  gives  a  comparison  ftl  cue  wire  current  with  that  in  the  free 
stream.  From  the  I  vs  vU      ..tion  data  on  tha  wtrt,  the  average 
velocity  as  a  function  o£  y  can  be  picked  off. 

Similarly,  a  measurement  of  wire  currents  at  a  position  first  with 
the  plate  heatec,  than  unhe-ated,  allows  calculation  of  the  temperature  rise 
above  ambient  by  the  method  of  Appendix  K. 


9)    »*«••; 


,  ja,*  .'.  I  -.•    6.   j       y     "  •:  ■ 


asi-v  s?<?   ri J.tv    bataoaiv  *.  UnwpseA. 


I 


* 


■ 


- 


b 

• 


1    *>*V3 


III.  RESULTS 

2 

A.  Figure  1  shows  uv/l?o»  as  »  function  of  vertical  distance  from 

the  plate  for  the  heated  and  ottheated  conditions.  Reasonable  agreement 
with  published  data  is  demonstrated.  The  value  of  this  quantity  is 
greater  in  the  heated  plate  condition  than  for  the  cold  plate,  as  was 
noted  in  reference  (1). 

Figure  2  is  a  graph  of  the  temperature   (Tg-Te)/(Ts-Te)  and 
velocity  (r/U»)  profiles  across  the  boundary  layer. 

Figure  3  is  a  plot  of  uv/U^,   *nd  vtfllao  <Ts-Te)  for  the  heated  plate. 
The  best  data  obtained  is  shown.  The  extreme  scatter  of  the  data  precludes 
confuting  fea  and  fih  as  functions  of  y  in  order  to  establish  an  analogy 
between  then. 

B.  The  basic  Instrumentation  selected  includes:  (a)  a  hot  wire  probe 
with  0.00015'*  tungsten  wires;  (b)  a   Flow  Corporation,  Hot  Wire  Anemometer , 
model  HVJBj  and  (c)  a  modified  Ballancine  Laboratory,  True  RMS  Voltmeter, 
model  320.  This  equipment  will  measure  accurately  the  properties  of  a 
momentum  boundary  layer.   It  is  probable  that  the  additional  properties  of 
a  thermal  boundary  layer  can  also  be  detected  with  this  instrumentation; 
however,  this  cannot  be  stated  conclusively.  Appropriate  recommendations 
are  included  in  section  VI.  The  equipment  is  illustrated  in  figure  4. 
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Figure   1 
Velocity  Correlations 
For  Heated  And  Unheated 
Plate 
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(a).   Equipment  Arrangement 


(b).   Probes 


Figure  4  -  Experimental  Equipment 
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IV.  &ISCU33I0H  ^F  RF,STOtS 

A.  The  plot  of  uv/U^e,  in  Pig.  1  lias  about  the  aaae  order  of  magnitude, 
although  somewhat  higher,  as  that  ahown  by  Johnson  in  ref.  1.  The  slope 
of  this  curve  is  such  greater  than  Johnson's  however.  He  did  obtain  two 
or  three  points  close  to  the  heated  plate  that  are  ssieh  higher  then  those 
for  the  cold  plate.  Further  irora  the  plate  Johnson* s  velocity  correlation 
quantities  are  almost  Identical  Lor  the  heated  aad  unheated  conditions. 

On  the  other  hand,  the  quantities  shown  in  Fig.  1  for  the  heated 
plate  remain  generally  higher  than  for  the  cold  plate  over  the  span 
■ensured,  this  difference  could  easily  be  due  to  experimental  error, 
however  there  is  a  possible  explanation.  Thermal  convection  from  the 
heated  plate  could  be  expected  to  disrupt  the  laminar  sublayer  causing 
some  increase  in  the  turbulence  of  the  flow  close  to  the  plate.  Fluid 
moving  from  the  sublayer  veulc  have  lower  velocities  and  therefore  tend 
to  correlate  a(~u)  with  a(-v)  to  sktid   to  the  quantity  (-uv)  for  the  cold 
plate. 

B.  The  plots  of  0/8 oo  and  (Tg-Te)/(Te-fe)  &a   a  function  of  distance  from 
the  plate  show  the  typical  characteristics  of  these  functions.  The  velocity 
profile  indicates  that  measurements  should  be  taken  closer  to  the  plate  to 
further  delineate  the  shape  of  the  curve.  To  do  this  a  much  shorter  hot- 
wire would  be  needed. 
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Extrapolating  the  curve  of  Tg-Te  in  rig.  10  to  the  plate  gave  an 
estimated  plate  surface  temperature,  Ts,  of  8  0  above  aubient.   Such  a 
low  temperature  difference  made  temperature  measurements  difficult. 

Another  more  accurate  method  of  measuring  temperature  would  be  to 
record  the  bridge  null  at  "cold  current"  (about  one  and  a  half  milliamps) 
while  traversing  the  boundary  layer.  For  this  a  bridge  galvanometer 
amplifier  is  required  to  enable  more  accurate  bridge  null  settings  than 
is  possible  with  this  model  &WB.  The  later  Wtk   models  do  Incorporate  such 
an  amplifier. 

C.  As  stated  in  the  Results,  the  scatter  in  the  computed  points  vt/Uoo  (Ts-Te) 
in  Pig*  3  did  not  warrant  calculation  of  can  and  Eh  by  equations  (3)  and  (4). 
This  is  one  of  the  most  important  ultimate  goals  of  these  measurements  that 
was  not  realized  in  this  thesis.  No  reliable  curve  of  vt  can  be  inferred  » 
and  several  attempts  to  deduce  better  values  of  this  quantity  produced  no 
better  results  than  those  shown.  A  higher  plate  temperature  would  likely 
give  larger,  taore-ea«ily-mfeaaure<l  quantities. 

The  data  does  not  sIjow  tiua  predicted  behavior  of  uv  and  vt  going  to 
zero  at  the  plate.  At  the  closest  approach,  the  center  of  the  wire  is 
still  0.068  inches  from  the  plate.  Any  decrease  in  these  quantities  must 
therefore  be  inside  this  distance. 

D.  The  equations  (69)  and  (70)  in  Appendix  B  for  two  wires  apply  also  for 
wires  mounted  on  an  X-probe.  Johnson  used  an  X -probe  with  both  wires  hot 
for  detecting  velocity  fluctuations  and  a  single  cold  wire  close  by  for 
temperature  variations.  The  hot-wire  theory  predicts  that  it  is  unnecessary 
to  use  more  than  one  hot-wire.  Since  equation  (71)  has  five  unknowns,  it  is 
necessary  only  to  make  measurements  with  one  wire  at  five  resistance  ratios 
and  solve  the  resulting  equations  simultaneously. 
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If  two  wires  on  an  X-probe  are  controlled  Independently  by  two  HWB*s 
with  outputs  measured  on  separate  meters,  equations  (69)  and  (70)  may  be 
solved  together  at  five  resistance  ratios  to  yield  the  boundary  layer 
properties.  The  two  Important  correlated  quantities  are  dominant  in 
equation  (72). 

The  practice  in  this  experiment  of  orienting  the  probe  at  -4 5°  simu- 
lated  an  X-probe  and  eliminated  the  problem  of  matching  the  two  wires  to 
close  tolerances.  Experience  with  the  V-probe  supports  the  comment  made 
by  Johnaon  (p.  54  of  ref.  1)  that  it  I*  unlikely  I1,  I,  <<0,  the  time 
constant,  lengta,  diataete*  and  c£u«  characteristics  of  two  wires  can  ewer 
be  siatched  simultaneously. 

A  few  comments  on  tae  V-pxobe  may  rfara  others  of  the  pitfalls  of  its 
peculiarities. 

It  was  first  discovered  that  the  resistance  of  the  common  current 
lead  and  center  needle  is  unavoidably  in  one  arm  of  the  bridge  circuit. 
Since  the  voltage  across  the  center  lead  ia  proportional  to  the  sum  of 
the  two  wire  currents,  the  bridge  aull  for  one  wire  is  affected  by  current 
in  the  other. 

The  second  problem  is  serious,  particularly  in  a  thin  boundary  layer* 
The  two  wires  are  sufficiently  separated  to  be  is  different  stream  condi- 
tions and  do  not  correlate  quantities  at  a  point. 

The  third  problem  has  been  discussed  -  that  of  matching  any  two  wires. 

Finally,  several  times  it  was  noted  that  wires  on  the  '/-probe  were 
not  taut  and  straight  several  hours  after  rewiring.  This  means  that  wires 
are  not  at  the  specified  angle  to  the  stream. 

B.   One  problem  that  detracts  from  the  accuracy  of  the  results  is  ambient 
temperature  variations  while  the  measurements  are  being  made  in  a  traverse 
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through  the  boundary  layer.  Using  the  minimum  number  of  two  resistance 
ratios  at  each  y,  it  required  four  to  six  hours  to  take  data  at  seven  points. 
Since  the  ambient  temperature  can  change  by  several  degrees  during  this  time, 
it  is  desirable  to  have  a  constant-temperature  air  input  to  the  tunnel. 

It  was  noticed  that  there  was  considerable  velocity  fluctuation  in  the 
free  stream  on  the  axis  of  the  tunnel.  This  masjfad  the  "e<Jge  ;  of  the 
momentum  boundary  layer. 

To  reduce  the  above  defecta  of  this  tunnel  it  is  desirable  that  the 
entrance  be  fitted  with  a  plenum  where  the  air  can  be  heated  a  few  degrees 
above  the  ambient  at  the  same  time  the  velocity  fluctuations  are  being 
esjeothed. 

Along  with  tunnel  entrance  modifications  other  test  facility  improve- 
ment? should  be  incorporated. 

Although  a  complete  survey  of  the  boundary  layer  should  include 
measurements  at  various  distances  along  the  plate,  it  would  be  well  for 
the  uninitiated  to  practice  in  boundary  layers  at  least  two  inches  thick. 
This  means  having  an  traheated  starting  length  of  at  least  40  inches  long 
and  a  heated  section  30  inches  long.  The  tunnel  should  be  three  to  four 
feet  wide  to  approximate  an  infinitely  wide  plate. 

The  heated  plate  should  be  well-insulated  on  the  bottom  and  equipped 
with  heating  colls  and  thermocouples  for  obtaining  3  constant  (q/A)  or 
constant  temperature  plate  as  desired. 

P.   Of  all  the  problems  encountered  in  this  experiment,  that  of  measuring 
accurate  (i.e.  lesa  than  51)  true  root -mean-square  voltages  caused  the  moat 
consternation.  The  best  available  information  indicated  that  the  Ball an tine 
Model  320  and  Flow  Corporation  TBM  were  the  only  two  meters  commercially 
available  reputed  to  measure  RMS  voltages  on  a   random  signal. 
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"hen  the  short  time  constant  of  the  Model  J2G  was  lengthened,  it 
became  the  most  reliable  instrument,  Uovever,  both  meters  exhibited  the 
annoying  characteristic  of  indicating  different  values  of  the  same  voltage 
on  adjacent  scales.  The  TBM  became  unusable  for  this  reason,  and  the  Model 
320  gave  ambiguous  reeding*  only  at  times.  Tike  limited  'crest  factor'  of 
the  Model  320  could  permit  peak  clipping  of  the  signal  on  the  lover  scale 
and  thereby  alter  the  «*eveforra.  However,  such  should  not  be  the  case  for 
the  T3M  with  its  much  higher  peak  factor,  " 

This  problem  was  not  resolved  and  is  the  most  likely  cause  of  the 
viste  scatter  in  the  calculated  results.  Knowing  that  this  problem  existed 
provided  the  hint  to  smooth  the  data  as  a  function  of  y  before  re-calculating 
uv  and  vt.  Even  so  the  improvement  was  not  marked. 

It  is  interesting  to  note  that  the  literature  reviewed  for  this  thesis 
makes  no  mention  of  true  RMS  measurements  or  time-averaged  readings  in 
similar  work. 

The  Model  320  could  be  improved  by  making  a  well-designed  permanent 
modification  to  obtain  a  long  time  constant,  mean  square  averaging  circuit. 
If  the  meter  time  constant  cannot  be  suitably  controlled,  then  a  five  to 
ten  minute  time  recording  of  the  output  could  be  graphically  integrated  to 
establish  the  best  value  of  each  voltage. 

A  more  elaborate  scheme  of  calibrating  the  meter  is  needed.  The  most 
appropriate  method  would  be  to  use  a  variable -frequency  generator  with  a 
calibrated  output  adjustable  from  one  millivolt  to  one  volt.  It  could  be 
used  to  resolve  the  problem  of  peak  clipping  of  the  signals. 

C.   One  teat  of  accurate  data  is  its  repeatability.  Time  permitted  little 
experimenting  along  this  line,  but  it  was  often  noted  that  the  bridge  could 
be  rebalanced  and  currents  reset  to  a  maximum  variation  of  one  part  in  150 
or  0.7X. 
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In  the  presence  of  signal  fluctuations  the  galvanometer  fluctuated  also. 
This  was  particularly  noticeable  at  high  resistance  ratios.   In  setting 
average  currents  it  became  necessary  to  average  the  needle  fluctuations  by 
eye.  Occasionally  this  led  to  variations  of  one  part  in  75  or  1.3%.  Pro- 
vision of  a  selectable  time  constant  on  the  galvanometer  would  assist  in 
setting  more  accurate  currents  enu  resistance  ratios. 

At  the  time  they  were  taken ,  meter  readings  could  be  averaged  within 
2%  in  spite  of  the  slow  fluctuations.  However,  it  was  rarely  possible  to 
repeat  the  data  at  a  point  because  ambient  conditions  had  changed.  On  the 
few  occasions  when  it  waa  possible  to  check  results,  meter  readings  were 
well  within  5%   herring  anomalies  already  discussed. 
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*•  CCWCLUSIOMS 

Hot  wire  anemometer  theory  predicts  that  both  of  the  desired 
correlation  quantities,  uv  and  vt»  can  be  obtained  with  a  single, 
matched,  X-wire  probe.  Due  to  the  extremely  small  volume  affected  by 
the  fluctuation  quantities  at  any  instant,  velocity-temperature  sensing 
elements  in  a  single  probe  is  preferable  to  separate  probe  sensing 
elements. 

The  instrumentation  system,  modified  as  recommended  in  section  VI, 
should  yield  data  of  sufficient  accuracy  for  calculation  of  all  the 
thermal  boundary  layer  properties.  The  results  presented  do  not  affirm 
this  statement  without  question,  but  they  do  indicate  that  further 
careful  Investigation  is  warranted.  Refinement  of  instrument  sensitivity 
and  experimental  technique  are  the  paramount  requisites. 
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VI.  RECCtttEWDATIONS 

The  following  modi f lest ion*  of  instrumentation  and  facilities  arc 
required  for  accurate  measurements  of  the  thermal  boundary  layer 
proper tia*: 

(1)  Permanent  Modification  of  the  meter  circuit  of  a  Ballantine  320 
to  obtain  a  long  time  constant  averaging  period; 

(2)  Accurate  calibration  of  all  the  range  scales  on  the  voltmeter  by 
use  of  an  accepted  standard  frequency  meter  which  has  a  calibrated 
voltage  output; 

(3)  Construction  of  an  X-wlre  probe  with  the  wires  oriented  accurately 

o 
at  90  to  each  other  and  their  resistance  matched  with  2X; 

(4)  Construction  of  a  longer  heated  section  for  the  plate  to  increase 
the  boundary  layer  thickness; 

(5)  Construction  of  a  larger  scale  wind  tunnel  with  an  inlet  heated 
plenum  chamber  to  insure  a  constant  temperature,  constant  velocity 
air  stream  in  the  tunnel; 

(6)  Procurement  of  an  additional  Wmt   modified  to  Include  a  sum  and 
difference  unit  in  lieu  of  the  amplifier  section. 

In  addition,  the  data  processing  techniques  should  be  refined 
as  follows: 

(7)  An  averaging  technique  must  be  used  over  and  above  that  which  is 
provided  in  (1)  above  -  e.g.,  averaging  of  sequentially  timed  ob- 
servations, or  graphical  integration  of  the  recorded  output  of 
the  voltmeter; 

(3)  All  five  resistance  ratios  permitted  by  the  wire  sine  must  be  used 
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to  calculate  all  fWe  tanas  of  aquation  #7/  until  tbair  order  af 
magnitude  ia  f  truly  established.  Programing  of  these  aquations 
for  computer  use  would  considerably  speed  data  analysis. 


APPEHDIX  A 


A  <  rose -«ec Clonal  area. 

c  Meat  e opacity  el  wire. 

C.  Meat  capacity  of  fluid  at  constant  pressure. 

C|tC2   Constant*. 

P  Wins  diameter. 

p  wire  diameter,  units  of  0.001  inch. 

a  Voltage  (general). 

«c     Output  voltage  from  the  amplifier  doe  to  se>sare»»ave  current 
alone  imposed  on  the  wire. 

ecn    Output  voltage  from  the  amplifier  due  to  square-wave  current 
plus  system  noise* 

e,.     Output  voltage  from  the  amplifier  due  to  system  noise. 

c9     Output  voltage  from  the  amplifier  due  to  velocity  fluctuations 
on  the  hot  wire. 

Output  voltage  from  the  amplifier  due  to  velocity  fluctuations 
plus  square-wave  current. 


aa  age  pies  system  noise. 

egg    Same  es  e«  plus  system  noise. 

ew     hot-wire  voltage  before  amplification. 

F      function  expressing  cooling  effect  on  hot  wire  by  velocity  and 
pressure  of  stream. 

7'     Defined  ••  rj   ,  the  derivative  of  F  with  respect  to  velocity  at 
constant  r.   p1*5 

f      Time-average  value  of  ¥. 
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I  instantaneous  deviation  value  of  F  from  F. 

I  Roto  of  boat  loos  from  tho  hot  wire  pox  unit  length. 

h  Fluid  enthalpy  per  unit  mess. 

I  Hot-wire  UC  current. 

I  Tlme-avorage  of  1. 

II  Defined  aa  rrr  1  ,  the  derivative  of  hot-wire  current  with  respect 
to  velocity  at  constant  pressure. 

i  instantaneous  deviation  value  of  1  from  1. 

ig  Instantaneous  deviation  value  of  hot-wire  current  due  to  Imposed 
square* wave. 

1^,1^  Hot-wire  BC  current  of  wires  a  end  b. 

Jj  Amplifier  multiplier  coefficient  of  ew. 

J2  *«,Ufi«  «,lttPU«  cffici-t  ot  *M  d.rivatl«  of  V 

K^  Meter  constant. 

k  Coefficient  of  thermal  conductivity. 

L  Length  of  hot-wire. 

n  Meter  reeding  of  true  RMS  voltage. 

*^»<W»Nn»l*s*Msc*inecn»|l(*n    Meter  readings  corresponding  to  voltages, 

e,  with  seme  subscripts. 

Mpy  Prandtl  number,  equals  cpy*7k. 

%#  Reynolds  number ,  equals  /±W/*>  . 

F  Fressure  of  fluid  in  the  boundary  layer. 

F  lime-average  of  P. 

p  Instantaneous  deviation  value  of  P  from  P. 

P*  Pressure  of  fluid,  units  of  atmospheres. 

q  Rate  of  beet  transfer  per  unit  time. 

R  Hot-wire  resistance. 

R  time-average  of  R. 
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r  Instantaneous  deviation  value  of  hot -wire  resistance  from  R. 

R9  Resistance  of  wire  in  equilibrium  with  ambient  fluid  temperature ,  T#, 

R€  Time-average  value  of  Re. 

r#  Inetantanaous  deviation  value  of  &#  from  R#. 

R  Resistance  of  wire  in  equilibrium  with  fluid  terapesaturc,  T  . 

*ft  «—««.  v.lw,  of  *g. 

r.  Instantaneous  deviation  value  of  R  from  Rtt. 
S  8      » 

RQ  Resistance  of  wire  at  reference  temperature ,  TQ. 

f  Temperature. 

T  Time-average  of  T. 

t  Instantaneous  deviation  value  of  I  from  T. 

Tm  Ambient  temperature  of  the  stream  entering  the  tunnel. 


T  Temperature  of  fluid  at  a  point  in  the  boundary  layer. 

Tg  Time-average  value  of  T  . 

T0  An  arbitrary  reference  temperature  for  computing  resistance  of 
wire  material,  handbook  value. 

Tf  Surface  temperature  of  the  heated  plate. 

Xw  riot-wire  temperature. 

X  Temperature  of  the  fluid  in  the  free  stream. 

u  Velocity  of  the  fluid  in  the  at-direction. 

mm 

»  Time-average  value  of  U. 

o  Instantaneous  deviation  value  of  U  from  U« 

u  Velocity  of  the  fluid  in  the  free  stream. 

V  Velocity  of  the  fluid  in  the  y-direction. 

V  Time-average  of  V. 

v  Instantaneous  deviation  value  of  V  from  V. 
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va»vb  Velocity  of  fluid  perpendicular  to  vires  A  and  8. 

V0  Velocity  of  fluid  at  angle    fron  perpendicular  to  wire. 

V  General  velocity  perpendicular  to  wire. 

v  General  deviation  value  of  Vp  fron  its  time-average. 

f*  Velocity  of  fluid  in  units  of  100  ft/sec, 

v  Velocity  of  fluid  in  ^-direction. 

V  Time-average  value  of  w. 

WA  Rate  at  which  work  is  done  en  a  particle  of  fluid  due  to  external 
fields. 

«g  Rate  at  which  work  is  dene  agsinst  viscous  stresses. 

X,  xt  Z    Body  force*  in  the  x-,  y-,  md   x-directione. 

x,  y,  a    distances  in  the  direction  of  mean  flow,  perpendicular  to  the 
plate  *&&  transverse  to  the  plate. 

**♦»  Temperature  coefficient  of  resistivity  at  T€. 

o(  temperature  coefficient  of  resistivity  at  f,. 
-•a  6 

o(Q  Temperature  coefficient  of  resistivity  at  T0. 

c<  Thermal  dlffusivlty,  equals  k/  f  c   , 

%  Resistance  ratio,  resistance  of  hot  wire  to  resistance  of  cold  wire. 

£k  coefficient  of  eddy  heat  (or  energy)  transfer. 

£n  Coefficient  of  eddy  momentum  transfer. 

&  Angle  between  perpendicular  to  wire  and  mean  velocity  vector. 

/*  Viscosity. 

tJ  Kinematic  viscosity. 

f>  Density. 

Z  Time. 
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Appendix  B 
BOUHDARY  LAYER  EQUATIONS 


The  general  equations  for  viscous,  constant-property  fluid  flow  ore 
quoted  below.  The  Navier-Stokes  equations  of  momentum  end  continuity  in 
rectangular  coordinates  are: 


(8) 


r  DX     da  * 


Continuity: 

&+frH+l7(pv)  +  li(ew)  =  0 


(9) 


The  energy  equation  is: 

where  w^  i§  the  rate  at  which  work  is  done  on  a  particle  of  fluid  by  an 
external  field  and  W|  ia  the  rate  at  which  work  is  dene  against  viscous 
stresses.  Also,  D/D?  represents  the  substantial  derivative  defined  as 

Ot    St   Zk       ay  T  ^T 

The  Kavier-Stokea  equationa  are  simplified  by  assuming  t wo- dimension- 
al ,  steady*  incompressible  flow  in  the  absence  of  body  forces  for  an  ideal 
fluid  of  constant  viscosity  and  thermal  conductivity*  Prandtl's  order-of- 
magnitude  treatment  further  reduces  them  until  the  result  is: 

e><     3y     f  3x      a^* 


< 
m»»  mmu$  mum  *m   ,3»*i»«f  s>* 


(  V-V)  J^l*  U*V 


^+X^3|-  =  ^, 
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[         si   ^  J  ~        1       ^6        xa  ^ 


°«(wi)£*W4  ♦(*!)■£♦* 


J&  +  JO*. 
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and 

iLL  +  _£V_  ,  o  C12> 

£>X     dv 

for  application  to  th«  boundary  layer  of  a  flat  plate. 

Similarly,  the  energy  equation  can  be  written  In  a  simpler  form  at 

ax         dy       pcP    a^«-       /pcT  \$x  / 

The  last  term,  m__  (  310  *  is  usually  considered  negligible. 
;pcP  Vdx/ 
Thue  the  three  equations  of  importance  are: 

Momentum: 

dx  ay  P   ax  ^l 

Continuity: 

Ml  +    -SV_  «  o  (15) 

6x  2>5 

Energy: 

i)ar_  +  v 'jl  =<£  a!I_  u*> 

with  boundary  conditions 


(a)    at    v-o   :      TaTs      or  M. %  <W     LUo 

3y  M 


(17) 


(b)   a*    y  =  oo  :     T=7oo      a/?i     lU  Uaa 

Multiplying  the  continuity  equation  by  U»  adding  the  result  to  the 
momentum  equation;  and  multiplying  the  continuity  equation  by  T,  adding 
the  result  to  the  energy  equation,  then  simplifying  both  gives 

<3x  *  3y   y     dy        I        / 

ax  ^  v     57  / 


o  .  _V£  +  JJ6. 

■  ■  .  >.     Ml 
Yj&)  ~-*U  +    J!4  _jU    IV.4  J£(J 


j^ci  y    ^T- 


mtam  \ 
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In  turbulent  flaw  Reynolds  assumed  thet  ell  quantities  ere  made  up 
of  en  average  velue  plue  e  superposed  fluctuation  component  such  thet 

U  =  U  +  u 

▼  »  V  +  V 

T  =»  f  +  t  (20) 

P-  •  t  P 
Substituting  these  la  equation  (18)   end  (19),  smiltlplying  out  end  taking 
tine  averages: 

_a_(?  +  pU*+f>i?)  *L(ac  30  -pUV-ptl^         (1) 

-£.  (Df  +  J*)  =  2_(<C  3L  -  Vf-Tl)  (2) 

In  general  uv  end  vt  nay  be  expected  to  have  values  other  than  aero 
es  ergued  below*.  The  U  velocity  profile  in  a  boundary  layer  over  a  flat 
plate  la  as  shown  in  the  sketch,  Fig,  5,  C hosing  a  point  at  x,y:   If  a 
•packet'  of  eir  goea  by  the  point  in  a  •+  V  direction,  elr  from  cloaer  to 
the  plate  will  be  forced  into  the  'hole*  left  by  the  packet.  On  the  aver- 
age, the  entering  eir  will  have  a  U-veiocity  leas  than  that  which  is  aver* 
age  for  the  point  x,y.  This  corresponds  to  a  (-u)  aa  observed  at  that 
point. 

Conversely,  movement  of  a  packet  in  the  (  -  v)  direction  gives  rise 
to  e  (+u)  velocity. 

Consequently,  where  there  is  a  gradient  of  If,  a  ( -+v)  produces  * 

(  -  u)  and  a  (  -  v)  produces  a  ( +  u) ,  therefore  the  product  la  negative  on 

the  average,  u  and  v  tarn   said  to  be  'cross -correlated*  such  thet 

■jlv  =   Limit    1   I     (llv)  at 

t-^oo  2.X   J  «  (23) 


as 
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Fig.  Si     Ve/ocJ+y    Profit* 


fig.  6-    Ve/oci+y  Correld+/on  Pro-f/'/s 
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Close  to  the  plate  the  fluid  becomes  laminar  in  the  sub-layer,  so  u 
and  v  tend  to  zero.  At  the  edge  of  the  boundary  layer  where  the  gradient 
of  U  disappears,  u  and  v  tend  to  cross •correlate  in  a  random  fashion  so 
the  average  is  sero. 

In  summary,  the  expected  behavior  of  uv  as  a  function  of  y  should  ap- 
pear as  shown  in  Fig.  6.  Parallel  reasoning  leads  one  to  expect  a  positive 
correlation  of  vt  where  there  is  a  gradient  of  T  such  that 

r-^o*  az  J_t  (24) 

These  results  are  sketched  in  Fig.  7. 

There  is  also  a  gradient  of  temperature  in  the  x-direction  producing 
a  ut  correlation,  but  usually  the  gradient  of  T  in  the  x-direction  is 
much  smaller  than  in  the  y-direction,  fluctuations  in  T  are  smaller,  so 
ut  is  expected  to  be  an  order  of  magnitude  smaller  than  vt.  Johnson  in 
reference  (1)  does  show  results,  however,  that  indicate  ut  to  be  about 
twice  the  magnitude  of  vt.  In  this  thesis  ut  cnacels  out  of  the  formulas 
used  to  calculate  vt.  Its  order  of  magnitude  compared  to  vt  Is  therefore 
not  determined.  t r 

It  must  be  remembered  that  near  the  'edge*  of  the  thermal  boundary 
layer  ut  and  vt  could  become  comparable;  however  the  methods  of  measure- 
ment are  unlikely  to  deteet  these  quantities  with  much  certainty  near  the 
free  stream. 

The  concepts  of  eddy  shear  stress  coefficient  and  eddy  energy  coef- 
ficient are  defined  by: 

*5  (3) 

and        vS  =  -  Ev,  &L 
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(4) 
where   £B  and  £^  are  probably  different   functions  of  y. 

When  the*©  definitions  are  substituted  in  equations  (1)   and  (2),   the 
result  is 


2  (T-hPU%p^)  =  2l 

57        ?      p  ** 


(25) 


and 


5k  djf  i         £y 


(26) 


Solutions  to  these  equations  have  been  obtained  for  various  simple 
geometries  and  assumed  empirical  distributions  of  pressure,  velocity  and 
temperature.  In  most  cases  the  ratio  of  £  m  to  £n  is  assumed  to  be  con- 
stant or  some  arbitrary  function  of  y.  Very  few  measurements  of  €m 
and  £h  are  available. 
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Appendix  C* 

HOT  WIRE  ANEMOMETER  THEORY 

A.  Steady  Plow. 

Hot-wire  theory  boo  boon  well  developed  for  various  Ideal  eoeee.  Moet 
of  the  theory  le  booed  on  King's  equation  which  appears  in  various  forms 
such  as 


I 


*<**  -  *e> 

or  alternatively 


=   c3  +  M  »Pr  hm'  U7) 


«» 
For  e  gas, p  is  proportional  to  ?/?»»  and  for  air  k  //**g  is  approximately 

independent  of  pressure  end  temperature  such  that 

H  =  (T„  -  Te)  F(FVp)  (29) 

where 


F(F»p)  s  [  0.58   f   6.6/p*?V  1 


(30) 


The  Flow  Corporation,  Cambridge,  Massachusetts  supplies  hot-wire  an- 
emometer equipment  which  nakes  use  of  several  interesting  techniques. 
First*  the  equipment  is  of  the  'constant  current'  type.  Second,  all 
measurements  are  wade  with  the  hot  wire  at  a  selected  resistance  ratio 
above  the  cold  wire  resistance.  Third,  measurements  of  unknown  signals 
from  the  wire  are  made  relative  to  a  known  square-wave  signal  applied  to 
the  wire  from  an  internal  circuit. 

A  summery  of  the  principles  of  wire  theory  will  give  adequate  under- 
*  This  appendix  is  bases  primarily  of  references  (3)  and  (12). 
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standing  of  wire  operation. 

For  a  hypothetical  fluctuating  flow,  the  power  Input  to  tha  wire  bal- 
aneaa  tha  power  carried  away  by  tha  stream, 

A  r  I  I,  =  (Tw  -  T#)  L  F(FVp)  (31) 

Alao,  neglecting  second  order  terns, 

*=R0[l    4  eC©(Tw  -  T0)J  02) 


wheTe  R0,  <£0,  and  TQ  are  arbitrary  corresponding  quantities.  These 

be  chosen  aa  R«,  <£e  and  T«  at  the  'cold*  temperature  of  the  wire  assumed 

to  be  in  equilibrium  with  the  streasi  ambient  temperature.  Therefore 

*  =  R#[l  +  <£e(Tw  -  T«)]  (33) 

Sowing  (33)  for  (Tv  -  T#)  and  making  use  of  ^=  K/Re»  equation 
(31)  becosjee 

X2      (1  -  -i-) U-  (34) 

r         cfe  *• 


In  steady  state,  therefore,  a  smasurement  of  I  will  give  the  9VP 
product  since  all  other  quantities  are  fixed.  In  practice  for  a  partic- 
ular wire  at  low  velocities ?n  calibration  of  I  vs  V  is  made  assuming  F 
is  nearly  constant.  Because  it  usually  results  in  a  straight  line,  the 
calibration  is  usually  plotted  as  I  vs  "]/Vp  giving  graph  that  may  be 
represented  by  an  .sanation  of  the  form 

I*   »  %  4  C2fT  (33) 

The  derivative  with  respect  to  velocity  Is 

«*    «       C*  (36) 


dF 


? 


rrF7 


In  Flow  Corporation  literature f the  quantity  dl/d?  U   demoted  by  I* 
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1.  Turbulent  Flow* 

When  the  stream  la  turbulent,  no  steady  state  exists  end  the  opere- 
tta* of  the  wire  must  he  analyzed  Iron  e  tins  dependent  approach.  The 
wire  stores  heat.  Since  it  cannot  change  heat  storage  Instantaneously, 
a  power  loss  balance  for  the  wire  is 

aTw  M<r~n 

From  equation  (33) 

$  T„  I       *  R 

— -   =•    — (3$) 

**  <£.  R«     R 

Take      t   -   ?.   +  r  (39) 

F(P?p)  =   F(PFp  +  Pvp)  (40) 


Substituting  these  in  conation  (37),  expanding  and  subtracting  the 


WP  +  %>  =  ?<*V  +  -r-  J  VF  =  F(*V  +  *  vp 


steady  state  equation  (31)  leaves 

c  A  I  -^—    =  «£«  R«  I2r  -  L  (R  -  R#)  rSp-lh  (42) 

Since  the  voltage  across  the  wire  is 

e^     -      1  r  (43) 

equation  (42) 
c  A  L 


j}y-'    =[^eRel2  -  "■»]  «w  -  L  <*  •  *e>  l  F*  *p  <**> 


Solving  lor  *'> 


L 


C  A  I*  d  fiv 


e     + 
w 


(R-R*)    IF*  w  1.  (R  •  R.)    IF*  d  r 


(45) 


Thus,  the  fluctuating  velocity  component,  vn,  is  e  function  of  the 

r 

voltage  and  the  tine  derivative  of  the  voltage  induced  in  the  wire  by  the 
velocity.  The  Flow  Corporation  8MB  performs  this  operetien  such  that  the 
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output  voltage  Is 


c  e. 


>«   =   Jl  «w  +  J2 


(46) 


r 


J2  is  a  constant  but  Jj  is  adjustable  so  that 


j   —   J*  (  «-»■-—«■»■«■-■»■■) 
1       *        c  A  L 
Solving  (45)  for  sw  or  d  e^V  df   and  substituting  In  (46) 


(47) 


r  jx  h  d  - 1^)  i  »•  i  r  j2  i  <r  -  r#)  p1 
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vn  (48) 


if  equation  (47)  is  true.  The  adjustment  of  the  ratio  of  J^  /  J«  it  made 
by  observing  the  amplifier  output  on  an  oscilloscope  when  a  square-wave 
current  is  fed  to  the  wire.  The  process  of  'compenseting  the  amplifier  is 
detailed  In  Appendix  B.  At  the  point  of  proper  compensation  with  the  wire 
in  still  air,  the  wire  response  to  the  square-wave  current  is 
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+      i. 


The  ratio  of  ett  to  ec  is 
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It  meet  be  realised  that  in  the  derivation  of  all  equations,  fluc- 
tuating i  usqiijnsiits  are  smell  compared  with  the  average  values.  Thin  al- 
Iowa  neglecting  second-order  products  of  the  fluctuating  component*  aa 
being  Immoniurably  small.  Even  -hen  this  la  net  strictly  true,  the 
•at hod  of  normalising  represented  by  equation  (51)  compensates  for  some 
deviation* 

C.  turbulent  Flow  with  Temperature  Fluctuations. 

A  fluid  in  turbulent  flow  past  a  heated  surface  will  become  heated 

in  none  fluctuating  manner.  The  effect  of  varying  fluid  temperature  may 

be  regarded  as  a  change  in  the  'cold*  resistance  of  the  wire.  £0,  such 

that 

*g  -   *g   +   *g  <*> 

=   *€  [  1  +  ie  <*g  -  *e>]  <55> 

Combining  these  with  the  power  balance  equations  in  the  presence  of 
velocity  fluctuations ^plus  the  relations 
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Noting  that     I   _   Tg  +   *  (57> 


(351 


and 


leads  to 


*g   =   *e[  *  +  <£,<Tg-Te)l 

*,   =   <£ e  *e  »  <58> 


2  i  +   2l'v i- -  t  s  0  (59) 

e8  and  ec  have  been  taken  as  the  outputs  from  an  ideal  amplifier. 
An  actual  amplifier  will  have  an  output  of  (e8  -f  a^)  and  ^c  >  eA  when 
this  output  is  fed  to  a  true  RMSvoltmeter,  the  latter  will  have  a  reading 

or    Men  =    «lVrJ^  1   <•«  +  V*  **"  <«> 

/*—  tt     J-t    1 Z__    -_     _ 

But  the  time  average  of  the  voltage  is  (et  <+-  en)/  =  ej  .+.   eg 

since  the  time  average  of  e^e^  —  0.  That  is,  there  is  no  correlation 

between  noise  and  signal. 

Therefore , 
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*\    +  4  (62) 


end  I         *» 


It  follows  that      y^n  ~     «£       =        ^y   .  (64) 

the  RMS  value  of  es.     For  the  square-wave 
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(65) 
Squaring  equation  (59)  and  making  use  of  the  meter  relations  gives 
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A  set  of  simultaneous  equations  can  b«  obtained  by  measuring  all 
quantities  at  several  resistance  ratios.  For  more  than  three  resistance 
ratios  the  unknowns  are  over determined  and  the  method  of  averages  or 
least  squares  can  be  used  to  evaluate  the  unknowns. 

B.  Resolution  of  Velocity  Components. 

The  foregoing  analysis  assumed  that  the  velocity  was  perpendicular 
to  the  axis  of  the  wire.  Within  70°  or  more  from  the  perpendicular,  the 
cooling  effect  on  the  wire  varies  as  the  cosine  of  the  angle,  that  is 
7g  s  v  cos  9  •  A  wire  at  45°  to  the  average  stream  direction  has  an 
average  effective  veloeity  1  f~\f~    times  the  stream  velocity.  This  ap- 
plies equally  to  the  steady  value  and  the  fluctuating  value  of  the  veloc- 
ity. 

Consider  only  the  fluctuating  components  in  the  turbulent  stream  y 
because  the  amplifier  filters  out  all  DC  wire  response  components.  For 
two  wires  at  ±  45°  to  the  stream, 

u  +  v 

▼a  =  —P=i («7) 
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Setting  va  s  v  ,  then  vfe  a  v_,  equation  (66)  becomes 
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Fortunately,  ch«  two  most   Important  quant I tie*  of  interest  la  the 
boundary  layer  analyst*  are  the  doalnant  teres  la  this  equation.   It  re- 
duces to 


*  I  «Jn  *  *2 
(2  i8> 


*su      "n 


i^scn  -  Msn/b  \Mscn  *  Hi 


=  —  8  I»2~5v 


(72) 


for  two  ideally  matched  wires  or  e  single  wire  turned  successively  to  the 
two  orientations  of  —  45°  to  the  stream. 
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The  experiments  were  conducted  in  •  tumicl  whose  test  area  has  a 
on*  foot  square  cross  section.  Atmospheric  air  was  drawn  through  the 
tunnel  by  an  exhaust  fan  powered  by  a  one  quarter  horsepower  AC  motor. 
rhe  flow  rate  tarougn  tne  test  section  was  controlled  by  a  damper  in  the 
exhaust  line.  The  maximum  velocity  attained  was  approximately  twenty-six 
feet  per  second,  lue  velocities  were  measured  by  a  pilot  tube  probe 
located  on  the  axis  of  the  tunnel  three  feet  upstream  from  the  test 
section.  A  micremonometer  was  connected  to  this  probe. 

The  boundary  layer  was  developed  over  a  flat  plate  which  consisted 
of  a  twenty-nine  inch  unites  ted  starting  length  and  a  twelve  inch  heated 
section.  The  heating  element  provided  approximately  constant  heat  transfer 
per  unit  area.  The  blunt  leading  edge  of  the  unheated  plate  tripped  the 

itum  boundary  layer.  As  shewn  in  Appendix  fc,  the  thickness  of  the 

ttum  and  thermal  boundary  layers  at  the  test  station  were  estimated 
to  be  1.32'  and  0.48'  respectively. 

(b)  Tjqhi  fjg  ftjUllvollmaftex, 

An  instrument  that  measures  the  true  root-mean-square  voltage  of  a 
random  input  signal  was  required.  There  are  two  commercial  meters  avail- 
able which  perform  this  function.  They  are  the  Ba I Ian tin*  Laboratories, 
Inc.,  model  330  and  Flow  Corporation,  model  TIM.  The  principal  difference 
between  them  is  the  manner  in  which  the  input  signal  is  averaged.  The 
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nodal  320  uses  en  amplifier  end  e  squaring  circuit  of  the  segmented  type. 
The  model  TIM  employs  a  unique  design  thermal  element.  The  principal 
characteristics  of  these  instruments  are: 


Characteristic  jsgjj  320 

meter  renge  (raw  volts)  0.0001-320  0.0001-250 

frequency  reapunse  (cycles  per  second)  15-130,000  2-250,000 

linearity  over  frequency  range         -  31  -  2X 

time  constant, including  meter  move-     2.23  15 

stents  (seconds) 

peak  factor  i  ratio  of  period  to      ■ "  ■  10 

pulse  width 

crest  factor  z     ratio  of  peas:  to  rms     3  

amplitude 

Due  to  the  presence  of  many  low  frequency  components  in  the  turbulence 
velocities  being  measured ,  a  relatively  long  averaging  period  is  required 
for  accurate  results.  The  problems  encountered  with  these  meters  are  pre- 
sented in  Section  IV. 

<c)  PT.obes 

Both  a  single  wire  and  a  v-wire  probe  mere  used  in  this  series  of 
experiments.  Each  had  tungsten  vires  of  0.00015"  diameter  with  sensitive 
lengths  of  0.040".  These  probes  are  shown  In  figure  4. 

The  single  wire  probe  was  used  to  simulate  an  X-wire  probe  as  noted 
in  the  discussion  of  results.  The  best  results  were  obtained  in  thie 
manner.  The  V-wire  probe  had  a  90°  included  angle,  and  the  sensitive 
segments  of  the  wires  were  located  as  close  as  possible  to  the  apex.  The 
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resistance  of  Che  vires  was  matched  to  less  than  two  percent.  The  probli 
previously  discussed  severely  Halted  the  use  of  this  probe. 

The  probes  were  mounted  on  a  traversing  mechanism  with  a  micrometer 
adjustment  to  accurately  measure  the  distance  from  the  heated  plate.  A 
calibrated  index  was  used  to  change  the  angular  orientation  of  the  probe. 

A  carefully  constructed  X-wire  probe  correctly  oriented  to  the  free 
stream  velocity  vector  would  satisfy  the  assumptions  made  In  the  derivation 
of  equation  (72).   It  would  eliminate  the  principal  problems  encountered 
with  the  V-wire  probe  -  i.e.,  apace  resolution  and  needle  resistance.  Such 
a  probe  was  not  constructed  for  this  thesis  due  to  fund  limitations. 


There  are  three  principal  circuits  required  to  transform  the  wire 
signal  into  the  data  required  by  equation  (72)  These  are  the  resistance 
bridge,  square-wave  circuit,  and  amplifier.  They  are  conveniently  pack- 
aged in  the  Flow  Corporation  £WB  unit,  fro  other  cooaer  dally  available 
anemometer  equipment  meets  the  requirements  for  measuring  fluctuation 
quantities. 

The  bridge  circuit  provides  a  method  for  accurately  measuring  the 
wire  resistance.  This  is  the  fundamental  quantity  which  is  sensitive  to 
velocity  and  temperature  fluctuations.  The  resistance  of  the  wire  at  amy 
eperetlng  condition  la  calculated  from  measurements  of  current  (I),  voltage 
<CW)  and  resistance  ratio  (  Jf  ) » 

The  amplifier  section  boosts  the  AC  components  of  the  voltage  across 
the  hot  wire  to  e  readable  level.  It  else  takes  the  time  derivative  of 
the  wire  voltage  and  adds  it  to  a  fraction  of  the  wire  voltage.  This 
fraction  ia  selected  by  adjusting  the  compensation  dial  to  obtain  minimum 
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signal  distortion.  To  properly  set  the  compensation,  •  square  wave  is 
impressed  across  the  bridge  circuit  and  the  amplified  wire  signal  is 
displayed  on  an  oscilloscope.  Correct  compensation  is  obtained  when  the 
square  wove  appears  linear  with  a  sharp  peak  at  the  leading  edge.  This 
vertical  peak  superimposed  on  a  square  wave  represents  the  addition  of 
the  derivative  to  the  input  signal. 

Any  random  input  signal  to  the  amplifier  nay  be  represented  as  a 
Fourier  series.  The  output  voltage  is  therefore  a  function  of  frequency. 
For  accurate  measurements  all  of  the  frequency  components  should  be  passed 
by  the  amplifier.  The  HW»  has  a  quoted  bandwidth  of  2  cps  to  100,000  eps. 
A  farther  error  would  be  introduced  by  overdriving  the  amplifier.  To 
prevent  this,  the  attenuator  dial  was  set  at  'is*  which  corresponds  to 
1/16  of  the  maximum  gain. 

The  square  wave  circuit  provides  a  reference  voltage  to  the  wire  for 
two  purposes.  First,  the  square  wave  is  employed  co  set  the  amplifier 
compensation.   In  addition,  it  is  used  to  calibrate  or  normalise  the  wire 
voltage  that  results  from  stream  fluctuations. 

One  other  piece  of  anemometer  equipment,  the  rfki,  was  used  to  s 
limited  extent.  It  contains  a  wire  current  control  circuit  and  e  eum- 
and-differenee  circuit.   It  is  used  in  conjunction  with  the  HUB  unit  and 
a  dual -wire  probe.  The  sum~snd-difference  of  the  instantaneous  signals 
from  two  wires  is  required  to  calculate  fluctuating  velocity  components. 

|t)  Oscilloscope 

The  output  signal  of  the  anemometer  amplifier  was  displayed  on  an 
oscilloscope.  A  Dumont,  type  322,  dual  channel  instrument  was  available 
in  the  laboratory.  The  duel  beam  presentation  will  be  found  particularly 
useful  with  an  x-wire  probe. 
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APPENDIX  E 
DATA  CALCULATIOM 

*•  Wjy  Calibration 

Before  each  wire  was  calibrated,  it  was  aged  at  a  lieatlng  current 
slightly  above  the  maxiwura  operating  range.   Its  velocity  response  was 
determined  by  calibration  in  the  free  stream  of  the  tunnel.  The  wire  was 
oriented  perpendicular  to  the  flow  and  measurements  were  made  of  wire 
current  for  various  resistance  ratios.  The  air  velocity  was  controlled 
by  a  damper  in  the  tunnel  exhaust  line* 

A  pitot  tube  with  a  micromanomater  was  used  to  determine  the  air 
velocity*  A  light  oil  of  0.327  specific  gravity  filled  the  manometer 
tube*  Ail  calculations  included  corrections  for  atmospheric  pressure 
and  temperature  in  the  tunnel;  however,  since  the  temperature  variations 
were  small,  no  correction  was  made  for  the  density  changes  of  the  oil* 
The  following  equations  were  employed  to  calculate  velocities  in  the 
tunnel: 

U  *    V^2gh     *     J2<g)(h)     0.827  /°M~|    * 

Pk    J 

.    I  2(32.2  ft/sec)     h  inchas  Q.827  x.  62.2  lbs/Ct3  I  * 

L  12  inches /ft  /°A  lbs /ft 3  J 


16.6 


W 


**•*•»    <°A  *     /°7triL<0»37g9)(  p  mm  Hg)  (0.07657)(0.3789)(p) 

<t°€  +-273.2*)  <*  *  273.2) 

T  (t  *  273.2) 
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Typical  calculation*  for  cwo  resistance  ratios  with  data  takan  as 
7/3 /S8  follow: 

pA  M    762.15  nH|i     tA  •     27.5°  C 

p.    n     <0,Q29plH762tl5)     .     ?J,UP       ,  0  07352  ibs/ft3 
'  A  (27.5  ♦  273.2)  300.7  °'0/3"  ™/« 

U     •     16.6  |    P-    ,,  •       ftfri6        \Th  *     61.2  i/T 

[o. 07352 J  0.2712 

h  fh  f"  I  (»a)  X2  <«a)2 

in  (in)*  ft/sec    fct/W    Jr«1.2     JT»X.4  V-  1.2       *«  1.4 

0.195  0.442  27.1  5.20  35.06  45.35  1229  2057 

0.153  C.391  23.9  4.39  34.43  44.66  1139  1996 

0.139  0.373  22.8  4.78  33.90  44.28  1149  1961 

0.117  0.342  20.9  4.57  33.80  43.75  1142  1914 

0.082  0.286  17.5  4.18  33.18  42.88  1101  1839 

0.040  0.200  12.2  3.49  31.35  40.75  983  1661 

0.000  0.000  0.0  0.00  23.25  30.32  541  919 

This  data  ia  plot tad  in  figure  8.  The  curves  follow  the  general  shape 
predicted  by  the  King  aquation  for  hot  wires: 

I2  «  Cx  ¥  c2  fo" 

2*  Se*er»ijftat;toa  oj  I'  =  9;pU 

The  alopa  of  the  curve  of  wire  currant  versus  incident  velocity  (I*) 
is  one  of  the  quantities  required  for  the  solution  of  equation  72.  Differ* 
entiatiug  the  King  equation,  we  obtain 
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where  ^2  *■■  til*  slope  of  the  calibration  curvet  obtained  in  part  1,  above. 
For  the  frae  stream  velocity  chosen  Cor  a  particular  experimental  run,  the 
value  of  (  \Tu)  Is  computed  and  the  corresponding  current  (I)  is  read  from 
the  calibration  curve.  With  thea»e  values,  the  required  slope  (I*)  is  cal- 
culated, keeults  are  shown  fot  a  tunnel  velocity  of  twenty-six  fe«t/second. 


ft/see     (ft/sec)*      •    (ma)Vft/sec)*     ma        ma/ft/sac 


26.0       5.10       1.2        132.6       34.86       0.1363 

1.4       217.4      45.03       0.2366 


3.  Amplifier  Cj 

At  dlaeussed  in  appendix  C,  a  square  wave  Is  employed  as  a  reference 
voltage  to  determine  the  anemometer's  response  to  e  fluctuating  voltage 
across  the  wire.  The  if^e  amplifier  compensation  is  set  by  observing  on  an 
oscilloscope  the  waveform  which  results  when  a  square  wave  is  impressed  on 
the  wire  circuit.  An  undistorted  uaveform  indicates  the  proper  calibration. 

To  determine  the  effect  of  compensation  setting  on  the  voltmeter  output, 
measurements  were  made  in  the  fre«  stream  with  the  hot  wire  oriented  verti- 
cally. The  values  of  Mgn  and  Macn  were  recorded  m   the  compensation  was 
varied. 
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64 
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m* 
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15850 
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15358 

0.0278 

This  data  is  plotted  in  figure  9.  The  response  is  approximately  linear 
with  no  eharp  will  or  peak;  hence,  proper  calibration  it  difficult.  It  ie 
estimated  that  the  error  introduced  in  the  cjientity  (Hg/Mj.)  when  setting 
the  compensation  is  r>t   the  order  of  8.2%  when  the  dial  is  varied  by  one 
dial  unit. 

Reference  12  gives  the  following  formula  for  computing  the  square 
current  of  the  tft«: 


3.oS 
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1220 

«?P 
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i  *    WiU 
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400  4  (w*)<  n   ; 

| 

where  (8*0  is  the  bridge  null  setting-  It  Is  important  to  maintain  the 
battery  supplying  this  current  near  „ts  fated  voltage. 
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5.  Temperature  Profile 

The  pertinent  formula  for  calculating  the  difference  <  A'O  bet 
the  temperature  (T„)  at  a  point  in  the  thermal  boundary  layer  and  the 
stream  ambient  temperature  is  derived  as  follows: 

*<p*>  *    JjLffcli<-JLA    *    *2  &p*q   f    *    ) 

2 

from  which,  (  h)      -   *  ,-,  *&   ■  *Z2jL. 
1  I2  /      R  -  R.      *  *  T , 

where  (Ij)  la  the  wire  current  for  a  cold  plate  and  (I2)  is  that 
current  for  the  heated  plate. 

Taking  R  *  Rg  1 1  +  <<    <1 ;  -  Tg)l  and  rearranging,  X  «  T g  +  (  ~^  ) 
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and  sub*  ti  tut  lag,   it  can  be  shown  that 


<*g  -  T«)     =  rf ; =f  (73) 


°|a-l>  [w^/i,)^  j 


For  an  anbient  temperature  of  24.8°  0*  the  coefficient  (  °^ •  )  Is 
computed  as  follows: 

c<  -   ?lfi .      y.tW»» .  0.004*4 

1  *«<0(ta-'i'o)        *  ♦  *>.00454(24.     2CTV 

Since  the  sire  i*  mora  sensitive  to  temperature  variations  at  low 
currants,  tha  data  for  the  temperature  ;>rofiic  was  taken  at  a  resistance 
ratio  of  1.1.  That  <lata  and  the  values  computed  %>ttb  equation  73  srs 
shown  in  the  table  below  and   plotted  in  figure  10. 
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The  norm* Used  temperature  profile  la  plotter  in  figure  2  of  faction  ill. 
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The  velocity  profile  across  the  momentum  boundary  layer  Is  obtained 
with  the  plate  cold  and  the  wire  oriented  vertically.  The  values  or  current 
recorded  are  then  used  as  arguments  to  enter  figure  b  and  determine  the 
velocity*  The  caea  and  results  are  shown  below  together  with  a  plot  of  the 
profile  in  figure  11.  The  normalised  plot  appears  in  the  *£SULTS. 
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7.  Boundary  Uiyer  ihicamqea 

The  following  equation  may  be  used  to  predict  the  momentum  boundary 
layer  fciiickne*6, 


£,.*  0.3?n(^-) 


Considering  a  feint  twenty-nine  inches  £r<xe  the  leading  edge  of  the  plate 
and  a  free  stream  velocity  of  tweuty*ix  feet  per  second,  the  thickness 
will  be, 
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manner  employing  Che  following  equation  tai&ea  iroca  reference  14, 
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8.  uv.  Correlation  Aa  "Cold"  Momentum  Boundary  Layer 

In  the  cold  plate  condition,  equation  7  nay  be  used  to  calculate  the 
quantity  uv. 


UV  "  fe 


ft-<      )        -( 


wL  -  H^   \     /  ML  -  M2 
sen"  **In 
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b 


Data  was  taken  at  two  resistance  ratios  and  averaged.  The  values  of  (I*) 
were  taken  from  section  (b)  of  this  appendix,  while  ig  was  calculated  for 
a  bridge  null  setting  of  483  by  the  method  of  section  4.  The  remaining 
data  is  shown  in  tabular  form  together  with  the  results  obtained,  which 
are  plotted  in  the  RESULTS. 
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••  SX  !S£  Xi  Correlations  in.  Thermal  Boundary  Layer 

On  the  hot  plate  condition,  equation  72  may  be  used  to  calculate 
the  quantities  uv  and  vt.  Rearranged, 
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**scn* 
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Data  was  taken  at  two  resistance  ratios  as  for  the  previous  calculations. 
Great  difficulty  was  experienced  in  obtaining  reasonable  results  because  of 
the  scattering  of  the  data.  To  reduce  the  random  errors  present  the  data 
was  plotted  in  figure  12  through  14,  and  faired  into  smooth  curves.  The 
remainder  of  the  calculations  were  carried  out  with  the  smoothed  data.  The 
raw  data,  the  faired  data  and  the  final  results  are  tabulated  below.  Normal- 
ised curves  of  the  results  appear  in  figures  i  and  3. 
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Raw  Data,  con. 


0    •  i.4 

y 

(in) 

4Ia     4Iw 
(ma)     (ma) 

(mv) 

(M8t»)a 
(mv) 

(mv) 

b  (^acn^a 

W 

0.068 

163.0    169.2 

7.5 

79.7 

108.0 

126.0 

144.0 

0.118 

163.8    171.4 

ti 

78.5 

102.5 

125.0 

143.5 

0,168 

166.3    173.3 

1* 

64.7 

83.1 

119.0 

129.0 

0.218 

167.7    174.3 

H 

56.0 

68.6 

112. 0 

124.5 

0.268 

168.5    174.3 

II 

49.8 

63.6 

108.0 

114.0 

0.368 

169.8    173.6 

n 

44.5 

56.0 

106.5 

111.0 

0.468 

170.0    173.7 

» 

41.0 

49.3 

104.0 

110.0 

Faired  Data 

Y  *   1.2 

y 

(in) 

*a 
(ma) 

(ma) 

(Me)2 
(mv)* 

(«a)| 
(mv)* 

(Kc>2 
(av)2 

0.068 

30.01 

31.93 

1690 

3520 

2350 

0.118 

30.65 

32.51 

1620 

2410 

»i 

0.168 

31.18 

32.89 

1130 

1600 

H 

0.218 

31.60 

33.06 

770 

1100 

m 
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31.88 

33.13 
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H 
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32.23 

33.19 
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» 
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O.dSI 

i»8tj 

t.es 

£.\ 

s.eax 

0.£dX 

880.0 

2.£*I 

o.esi 

t.sox 

e.8* 

H 

« 

8.£dX 

8X1.0 

O.SSX 

o.ei: 

X.£8 

t.*8 

H 

£.£U 

£,ddX 

88X.0 

Mil 

O.SII 

8.88 

o.de 

M 

C.MX 

S.Vdl 

8XS.0 

o.*xx 

0.801 

1 ,  M 

g.t* 

H 

e.Mx 

«.88X 

8d£.0 

O.XXX 

e.aox 

o.de 

l.H 

M 

•4 

88E.0 

O.OXX 

0.*>X 

e.t* 

O.I* 

II 

■  u 

O.OU 

88A.0 

ftj«a  b*iiA% 

£.X  «  ^ 

£(3JO 
£(v*) 

d^; 

*(*«> 

.1 

(MB) 

v 
0*) 

Ntt 

0£C£ 

owi 

£«.!£ 

xo, 

,0£ 

880.0 

Of*S 

0S8X 

xe.££ 

•», 

>oe 

8X1.0 

00*1 

oexx 

W.S£ 

II 

it 

88X.0 

■ 

OOXX 

ou 

do.ee 

Oft. 

,it 

8XS.0 

■ 

m 

0«       ' 

£!.££ 

88, 

,X£ 

88S.0 

oee 

0£C 

«.€£ 

a 

.££ 

88C.0 

II 

OM 

HI 

0S.££ 

**< 

Jt 

88*.0 

60 


Faired  Data,   con. 


»     1.4 


7  la  lb  <*«>*  (M.)g  (M>2 

<*>  <««>  <~>  jZ)J  ^2  ffi2 

9400 


0.068 

40.87 

42.30 

6320 

11650 

0.118 

41.13 

42.81 

6040 

10250 

0.168 

41.55 

43.22 

4190 

7380 

0.218 

41.89 

43.42 

3130 

5000 

0.268 

42.10 

43.50 

2470 

3840 

0.368 

42.41 

43.56 

1830 

2790 

0.468 

42.56 

43.58 

1600 

2400 

Results 

y 

w  2 
(ft/sec) * 

co 

vt 
(ft  C /sec) 

vt 

(in) 

»«oCVV 

0.068  +4.69  4 0.00694  +7,96  40.0383 

0.118  -2.74  -  0.00405  +1.74  +0.00837 

0.168  -3.35  -0.00496  t-0.421  40.00202 

0.218  -1.41  -0.00209  *0.633  40,00304 

0.268  -0.641  -0.00095  ♦O^?  *0.00354 

0.368  -0.205  -0.00030  40.686  40.00330 

0.468  *0.002  40.00000  fO.693  fO. 00333 
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